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The general theory of the motion of a single rigid body through an infinite incom- 
pressible fluid is well known, chiefly through the work of Thomson and Tait* and 
KiRCHHOPF,t and we are able to calculate numerically the results in the case of the 
sphere, the ellipsoid, and a large number of cylindrical surfaces. The theory of the 
motion of two or more bodies in a fluid has naturally not made the same progress, 
and we are unable to determine the form of the expressions involved for the general 
motion of any particular solids. So far as I am aware, the first attempt was made by 
Stokes, in a paper read before the Cambridge Philosophical Society in 1843, entitled 
" On some cases of Fluid Motion. "J In this paper, amongst other problems, he con- 
siders the case of* two spheres. He determines the instantaneous velocity potential 
for two concentric spheres and for two concentric cylinders with fluid between them, 
and finds that the effect of the fluid is to increase the inertia of the inner sphere by 

a mass ■=.\.— — - of the mass of the fluid displaced, and that of the inner cylinder 

(./ ~— ~ (v 

Q& _1_ /v2 

by a mass — — - of the mass displaced, a, b, being the radii of the spheres or cylinders. 

He also approximates to the cases where one sphere is moving in the presence of 
another in an infinite fluid ; and also in the presence of a plane, the method used 
being first to calculate the velocity potential for any motion of the points of the 
plane, and then suppose them actually animated with velocities equal and opposite 
to the normal velocities of the fluid motion at those points if the plane had been 
removed. He applies the same method also to the consideration of the motion of 
two spheres. In a note in the Report of the British Association at Oxford, 1847, 
he states the theorem given by me in § 4 relating to the image of a doublet whose 
axis passes through the centre, and mentions that this will easily serve to determine 
the motion. In 1863 Herr Bjerknes communicated a paper to the Society of 
Sciences at Christiana, on the motion of a sphere which changes its volume, and in 

# Nat. Phil., p. 264, new edition, p. 330. 

f BORCHARDT, Bd. 71. 
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which he approximates for the motion of two spheres. I have not been able to see this 
paper, nor some others which he presented to the same Society at some later periods ; 
but he has given an account of his researches in the ' Comptes Rendus/* together 
with some historical notices on the development of the theory. He does not seem, 
however, to have been acquainted with the important paper of Stokes above referred 
to.t In 1867 Thomson and T ait's ( Natural Philosophy ? appeared, containing general 
theorems on the motion of a sphere in a fluid bounded by an infinite plane, viz. : that 
a sphere moving perpendicularly to the plane moves as if repelled by it, whilst if it 
moves parallel to it it is attracted. In a paper on vortex motion in the same year 
(Edin. Trans., vol. xxv.), Thomson proved that a body or system of bodies passing on 
one side of a fixed obstacle move as if attracted or repelled by it, according as the 
translation is in the direction of the resultant impulse or opposite to it In the 
'Philosophical Magazine' for June, 1871, Professor Guthrie publishes some letters 
from Sir W. Thomson on the apparent attraction or repulsion between two spheres, 
one of which is vibrating in the line of centres. Results only are given, and he states 
that if the density of the free globe is less than that of the fluid, there is a " critical " 
distance beyond which it is attracted, and within which it is repelled. The problem 
of two small spheres is also considered by Kirchhoff in his ' Vorles, u. Math. Phys./ 
pp. 229, 248. In his later papers Bjerknes takes up the question of cc pulsations " as 
well as vibrations. Of solutions for other cases than spheres, Kirchhoff has con- 
sidered \ the case of two thin rigid rings, the axes of the rings being any closed 

* ' Comptes Rendus,' torn, lxxxiv., p. 1222, &c. 

f Not only Herr Bjerknes, but several writers on the Continent seem to be unacquainted with this 
paper of Stokes, and also with Green's papers. Kirchhoff, in his 'Vorlesungen iiber Mathematische 
Physik ' (second edition, p. 227), says that Dirichlet first treated the motion of a sphere in a fluid in the 
Monatsberiehte der Berl. Akad.' in 1852, and Clebsch that of the ellipsoid in 1856, in c Crelle,' Bd. 52. 
Bjerknes also repeats the same statement, and Clebsch in his paper regards Dirichlet as the first to 
solve for the sphere. In his paper Dirichlet says : " Wie es scheint, ist bis jetzt fiir keinen noch so 
einfachen Fall der Widerstand, den ein in einer ruhenden Fltissigkeit fortbewegter fester Korper von 
dieser erleidet, aus den seit Euler bekannten allgemeinen gleichungen der Hydrodynamik abgeleitet 
worden." The fact is that Green in a paper read before the Royal Society of Edinburgh in 1833, 
entitled " Researches on the Vibrations of Pendulums in Fluid Media " (Trans. Roy. Soc. Edin. ; also 
published in the Reprint of his papers, p. 313), and written without the knowledge of Poisson's paper of 
1831, " Sur les mouvements simultanes d'un pendule et de Pair environnant," treated of the motions of 
an ellipsoid moving parallel to one of its axes. He obtains the velocity potential as an elliptic integral 
for a motion parallel to an axis, which also of course contains implicitly that for the sphere. He shows 
that it is necessary to suppose the density of the body augmented by a quantity proportional to the 
density of the fluid. For the case of the spheroids moving in their equatorial planes or parallel to their 
axes he completely determines this quantity, whilst for the sphere he finds that it is one-half the mass of 
the fluid displaced. The first place in which I have been able to find the well known form of the velocity 
potential for a sphere is in Stokes' paper of 1843 before mentioned. He obtains it as a particular case of 
a more general problem, and refers to it as the " known " value for the sphere, The equations of the 
lines of flow were, I believe, first given by Dieiohlet, 
+ Borchardt, Bel. 71, 
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curves and the sections by planes perpendicular to the axis being small circles of 
constant radii, and he arrives at the result that their action on one another may be 
represented by supposing electric currents to flow round them ; and I have recently 
solved the problem of the motion of two cylinders in any manner with their axes 
always parallel. The velocity potentials for the motion of the two cylinders are 
found in general as definite integrals, which, when the cylinders move as a rigid 
body, are expressed in a simple finite form as elliptic functions of bipolar coordinates. 
The functions involved in the coefficients of the velocities in the expression for the 
energy have a close analogy with those for spheres arrived at in the following 
investigation. 

1. Our first aim will be to find the velocity potential for the motion of the fluid in 
which a sphere is fixed and in which a source of fluid exists. By the image of the 
source in general is meant that collocation of sources or sinks within the sphere which 
produces outside of it a fluid motion which in conjunction with the original source has 
no normal motion across the sphere : in other words, that " mass " of positive or 
negative sources which produces across the surface of the sphere a normal flow equal 
and opposite to that of the outside source. When this " image " is found, the way is 
theoretically clear to finding the velocity potential when two spheres are fixed in the 
fluid, and thence, by distributing over the surface of the spheres sources proportional 
to the normal motion of the surface at that point, to determine the velocity potential 
when the two spheres are moving in any manner. In the case of an electrical point 
the image is, as is well known, a negative point at the inverse point of the other. In 
the case of fluid motion the image is, as will be shown, a positive source at the inverse 
point, together with a negative line sink stretching from this point to the centre of 
the sphere. 

Fig. 1. 




2. Take O the centre of the sphere for origin and let the axis of % pass through the 
source S. Let the radius of the sphere be a, and the distance of S from the centre 
be 6. Then the velocity potential will clearly be symmetrical about O S. The 
velocity potential for the unit source at S can be expanded in, the series 

J- J- X .qo / J-^ 

7 — 2,1 
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which holds good for points where r<b, whence when r=a(<b) the flow into the 
sphere at any point (0) is 

Expand the potential due to the sources, &c, inside the spheres in a series of 
spherical harmonics 

Hence the flow out of the sphere, for points just outside, is 

1 

and this must be equal to the other, whence 



?i /&^ +1 



Y " = ~TTiU' p » andY o=° 



and 



n-\-l 



% a? n+l 



(pr) 



n+l 



n 



Consider 






^l__ /y»_>^_p /X<a 



the potential for a source // at a point on O S inside the sphere at a distance \ from 
the centre. Then 

rv ju/\ l x w+1 



a 9. 



Comparing this with the expression for V, we see that if we make X=- and 



, a 
u=;X source 






6 -* r ^+i - - • x 



o ^ 



a 1 1 r A rf\ 

b ^/r* — 2\r cos + X 3 a- J ^/r 2 — 2Xr cos + X 2 



a 3 



?*.6. ? V is the potential of a source at the distance y- from O whose magnitude is equal 



MR. w. m. hicks on the motion of two spheres in a FLUID. 459 

d , ^ eft 

to - of the source at S, together with a line sink extending from O to the distance ■—, 

1 
the line density of the sink being - X source at S. 

Performing the integration for V, we find finally that the whole velocity potential 
for a unit source at S is 

• -L. - *%*«** JL W .1* 

^=—^13+ V=: 



SP ^ s/t*- 2br cos -f 5 3 b s/r* - 2\r cos + X 3 

. 1 i \—r cos + \/r 2 —2\r cos # + A, 3 
a ° r(l — cos 6) 

where X= - 

o 

It is easy to verify this value for <f> by direct differentiation. 

If we apply the same method to find the velocity potential for the motion of fluid 
inside a sphere under the influence of a source inside, the integral becomes infinite 
unless the source is zero. The case is of course physically impossible since if fluid is 
generated within the sphere it must pass through the boundary. But if we also place 
an equal sink at any point within, the motion is then possible, and the expression 
becomes finite. S being the source let S' be its inverse point with reference to the 
sphere, and S" any point on the line S S' produced to infinity. Then the " image 77 of 

S is a source ~~- at S' and a line distribution of sinks of line density - from S' to 

infinity. Let S x be an equal sink, then its image and that of S will produce potentials 
with finite derivatives. In fact, the potential at P will be 



^—^ 



1 1 ,a 1 a 1 1. OS'-rcos0 + ST l-cos0, 

■ — - — j — *■ «—» — • — — loo* 



SP S^ ' b ST b, S\V a B OS'i-r cos 0, + S\Y 1- cos 



where #, Q 4 are the angles O P makes respectively with S, O S^ 

3. The expression found for the motion when there is a single source outside the 
sphere enables us to deduce the velocity potential for a single sphere moving through 
an infinite fluid. Taking the direction of motion as the axis of x 9 from which we will 
suppose 6 measured, we may arrange a surface distribution of sources proportional to 
cos 6dS and integrate over the surface of the sphere, or we may employ the simpler 
method used in a paper in the ' Quarterly Journal of Mathematics 7 for March, p. 128. 
The first gives us the velocity potential when the sphere moves by an integration 
which would be laborious. The other gives directly the potential, when the sphere is 
fixed and the fluid moves past it, by means of an easy differentiation. Putting a 
source at #=& and an equal sink at #= — &, let these move off to infinity, increasing 
indefinitely as they do so, yet so that the motion at a finite distance from the origin 
is finite* In the limit we clearly get the case of fluid flowing past the sphere. 
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We have to find the limit when b~ go :-=& of 

6^ 



^>= —/a 



1 



yV 3 — 2br cos (9 + P \/r^ 4- 2&r cos 9 4- & 3 



a / 1 _ 1_ 

1 , X— r cos + \/r 2 -~2\T cos + \ 3 1 + cos 6 

I S"\ /"V» „_„- . , -__ . : 

"" to A, + r cos 6 + v^^+ 2\r'cosT+ A?'l — cos 



C6 



When 6 is large and X small this is easily shown to be 



Hence the limit is 



0= — fd 2r cos 6'+ — ^ — +r+ • • 



I 7 / . ( ^ X 

0=— A( 2ar+~j 



If the velocity of the fluid at an infinite distance parallel to x is ^ towards the 
origin, then 

Also impressing on the whole system a velocity u, the sphere moves with velocity u in 
an infinite fluid, and the potential function is 

, _ a?ux dhi cos 9 

The well-known form of <f> in this case. 

4. If now two spheres A, B are present in the fluid, and we consider the series of 
images resulting from the first image in A, we see that they very rapidly become 
extremely complicated, e.g., the first image is a source and line sink ; the image of 
this in B consists of (l) a source and line sink, (2) the image of the first line sink or a 
line sink (segment of a circle), and an area source bounded by this last line sink and 
two straight lines from the centre. It is, therefore, hopeless in this way to find first 
the velocity potential for a source in the presence of the two spheres, and thence the 
potential for any motion of the spheres. But now suppose A fixed and B moving in 
any direction* If A were not present the velocity potential of B would be that due 
to a doublet at its centre, whose axis lies in the direction of the motion of B. The 
effect of the introduction of A will be to produce a series of images of this doublet, 
lying inside A and B. This method dispenses with the necessity of integrating over 
the spheres when we have found the velocity potential for the doublet. In the special 
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case where B is moving in the line joining the centres, the image becomes simplified 
and reduces to a single doublet. For let us find the image of a doublet whose axis 
passes through the centre of a sphere. 

The doublet is formed by allowing an equal source and sink P, P' to indefinitely 
approach one another, their magnitudes increasing indefinitely, yet so that /x.PP' is 
finite. Now let P, P' lie on the line through the centre of the sphere, and let Q, Q' 
be their inverse points ; moreover, let the limit of /x,.PP'=&. Then the image of P, P' 

consists of a source ~p at Q, a sink ^ at Q', and a line source (supposing P outside 
F, and therefore Q' outside Q) along Q Q' with line density -, also the quantity 

-•QQ', together with the sink at Q', is equal and opposite to the source at Q, and we 

may suppose it added to the sink at Q', when they become equal Now as P, P' 
approach to coincidence so do Q, Q', and the image of the doublet k at P becomes the 
doublet at Q, whose magnitude is the limit of 



aP 



Qp*^{^ — Op'pP' — OP 3> 

/ a \3 

i.e., one of opposite sign and magnitude ("•--) X that at P. The same result can 

easily be shown to follow from the analytical formula in § 2. 

The case where the doublet has its axis perpendicular to the line joining the centres 

has more analogy with the case of a source. The image here consists of a doublet of 

the same sign at the inverse point, with a trail of doublets of opposite sign extending 

to the centre. 

Fig. 2. 




Let, as before, P, F be equal source and sink, Q, Q' their inverse points with respect 
to the circle. 

Then at Q, Q' we have a source and sink of magnitude ~, and in the limit we have 
a doublet 



^QQ 7 _^ 



a 



OP VOP 



t 

a 
MDCCCLXXX. 3 O 



Also, if B, B/ be corresponding points on O Q, Q', we have a line density —--at 
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V> 



R and -f - at R'. Consequently when P, P' approach indefinitely so do R, R', and 

CO 

we get a line doublet along Q, whose line magnitude at any point R is the limit of 



i.e., proportional to the distance from the centre. 



Fig. 3, 




5. Supposing that the positions of all the images of the doublets and their mag- 
nitudes are known when the sphere A is moving along the line B A, we proceed to 
find an expression for the kinetic energy. Let p u be the distance of the n tl1 image in 
A from A, and <r n the distance of the n th image in B. Also let the magnitudes of the 
doublets there be [i m v n respectively. Let <f> be the velocity potential of the motion , 
and <j> m <f/ n the parts of <f> due to fi, n and v n . Then denoting the kinetic energy by T 

2T= - f[> \fdS=z - 27ra?u\X<l>] sin cos 0d6 

where [</>] is the value of <j> at any point (a.ff) on the sphere. Now ^=2^»+2<£'» and 
the part of T due to <j> n will be 



2T= 



J {& 2 + 2p tt aeos0 + /VT 



2p n a/j,y 



Now 



rf 



+1 (p + a,ft)iid/A__ 
^i{c<? + p* + 2 pap,} 

1 



rff+i 



ytfrd^tf, 



cZpJ ^y'tt 3 -f /0 3 •+ 2/Xty& 



When p=/>^ p w <a and the above becomes 



zrd(p+ a )(p*+ a *-P a )-(p^\ a )(p*+ a2 +P a )} 
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Similarly when p~cr u <r n >a and it becomes 



d 2a 
dp'Sp?'' 



4a 
3p» 



and 



2T= — f ir.u%™ p n -\-^ira?u%\ 



oo / Vn 



<T 



n 



Also [jl the original doublet = 
sphere A 



ahi 



and if M 1 be the mass of fluid displaced by the 



2T=pi 1 # 



1+S 



<x f^n 
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But from what has been shown before 



w- 
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;r 



cr 



S^w 
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Hence 



2T=|M^ l + 327(- 



« * • 



• # • • • \ -^ / 



By § 4 we have, if c is the distance between the centres, 



(X 



w 



n 



cr 



^n 



n 



Also 



Hence 



I \3 
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«■ 
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fin-l 
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Again 



PnPn—\ • • • • Pi 



l\3n 

a/ I (c—pn^) .... (e — p-^c 
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whence 
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Put /o»=w»+#, and choose x so as to make the constant term vanish. To find x 

we have 

<x 2 -f c 3 — & 3 






St 



c 



x-\-a?=0 



Fig. 4. 




Now let C { C 3 be the inverse points of the spheres, and O the middle point of 
Cj Og. Put Oj 02 ==: 2a. Inen 

OA=^\ i +a i -^-^—^-=r 1 Ba,j ( 2 ) 



C = r 1 + r s 

P 
Further, P being any point on the sphere A, denote the constant ratio ~~ by q x 

and let q % be the similar constant for the sphere B. Then 



9i 



ft= 



\ + r l —a_\ + r 1 

X—Tj + a a 

X + r Q — 6 \ + r„ 



a 



r x -X 



r 3 -\ 



& 



S * 



(4:) 



The equation to determine x now becomes 

x^-—2r l x+a?=0 
The roots of which are ^ifcX. 

Choosing the positive sign, the equation of differences becomes 



a 



2^ 
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2^ 



K+ #i— ~~ K- 



U n U n ^ — \ ^2 " j w n~T\^l~~ I w «— 1 ' 



Now a s =a5 1 a? 2 whence writing - for u n we get 






^(c— ■ #y 
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Here 



and 



C "~~ Xa I a ~f" J\j 



2 



C "~~ JO-i in """" A> 
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•2 



2 T l 



Xo r-, —X * ] 



% 



Whence writing ~ ~q 
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x x (c-x^ (r l +\)(r^ + \) 
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But 
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Whence 



2T=iM>» 1+3(1 -qr") 5 ^ 
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We shall denote, in what follows, 



(i-gr 3 ) 3 2r r 



9 



n 



9l -y» 



6. If the sphere B is also moving- along 1 the line AB, the kinetic energy of the 
fluid will be of the form 



£J X — - Q i-V-i ] Cvl 



l + 3Q(^j 



+iM 3 % 3 { 1 +3Q(g 2 .g)} +L%% 



It remains, then, to find the value of L. 

It is easily seen that L depends on the part of <f> belonging to the images of B s 
motion taken over the sphere A, together with that belonging to the images of A J s 
motion taken over the sphere B. Let now dashed letters apply to the images, &c, 
of the B system, then using the results in § 5, the part of L due to the integration 
over xjl is 

— —-girVriZii LL ^-*T 3"7T<% U-iZiQ i ~71 

But as before, remembering that now the original doublet is in B, 
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a 2 



and determining A by the condition that pi=~~ we shall find 

c 



ci/JLlUL 



an 



d 



ac — p',^ ■* 1— g 3 * 



L 1 =-2flrw 1 %Si — : — 1 — s 



==~16^ 1 %X%f^£ A 



C 1— g* 

» \8 



Similarly L, 2 = same quantity. 

Therefore, denoting by M' the mass of fluid contained in a sphere of radius unity 

L = — 47m 1 t%Q'(g) = — SWu^Q^q) 
where 

Q'(gr)=2i [rz^nj • • * w) 

Tables for Q and Q' are given at the end of the paper for equal spheres, and for the 
case of a=2&. 

7, When the sphere A is moving perpendicularly to B A, the original doublet is one 
perpendicular to the line B A, as also its images. Suppose A is moving along the axis 
of x, A B being the axis of %. Then the normal velocity at a point P on the sphere A 
is v sin cos x> ( a «0«x) being the polar coordinates of P; and the kinetic energy is given 
by 

[re [%n 

2T=— a 3 v [<£] sin 3 6 cos x dddx 



Let fx be the magnitude of a doublet at a point distant p from the centre of A; the 
part of <f) depending on this is 

fir sin 6 cos % 
{r 3 + /> 3 + 2pr cos 0}* 

and the part of 2T depending on this is 
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The integral of which is 
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Writing v and cr for fx 9 p for doublets outside the sphere A, we obtain 



4 t 47ra 3 z>u 
^tt^v and -~ o ~ 
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whence 



2T= — M^S^j "^+7; 
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Now any v at the distance cr produces an image in A consisting of a doublet v I - 



a 

or 



a 



CI* 



at a distance — , together with a line sink stretching from this to the centre, whose 

line magnitude is X distance from the centre. Hence the whole amount of the 

aa 



image is 



\<r/ a<r\<T J * \(7 



Now every fi except [i forms part of an image of some v 9 and of that v only* 
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(7) 



The S extending to the whole mass of images inside A. 

8, If A has also a motion along B A, together with one perpendicular to it, T has 
no term depending on u 9 v; for it is clear that if the sign of v is changed, then the 
kinetic energy must be the same as before. 

If B moves also perpendicularly to B A, T will have additional terms in v 2 2 9 and 
v l9 % The coefficient of v% will be analogous to that for v* 9 whilst that for v l9 v 2 , as 
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in the case of u l9 u 2 , consists of two parts, depending on the integration over the two 
spheres. As in the case of L this coefficient 1/ is 



L\v fi =—i7rv l l l \i/+v 



a 



a 



3 



%7TV 2 t 



-*©'] 



dashed, letters referring to the motion of B, /x, v referring to images within A and B 
respectively, and cr' <t denoting distances from the centres of A and B. This may be 
reduced as in the former case to 

the /x' being the images in A of B's motion perpendicular to A B, and v the images in 
B of A's motion perpendicular to A B. 



1/^0= — 4:7rv l v % 1 



V* 



v 



47rV 2 ja S 



o 



V 



My 



where v Q9 /x, are the original doublets at the centres of B, A, i.e. 



l\c 



arv-. 



P, 



0' 



/*o= — 



2 



whence 



/*' 



L^fM^fJ-j+fM^ 



V 



M/ 



(8) 



in which last the ratios *--, — do not contain v l9 v 2 . 

In the case of general motion of two spheres, each will have three components of 
velocity, u l9 v l9 w Y \ %, v% 9 w 2 ; and, in general, the expression for the kinetic energy 
will contain 21 terms. In the case in question we can easily see that the coefficients 
of 12 of these vanish. For consider the term in u l9 w 2 — suppose v l9 w l9 u% 9 v% all zero; 
the energy, from the symmetry of the motion, must clearly be unaltered if we reverse 
the direction of w 2 . And this can only happen if the coefficient of n^w^O. In this 
way we find the terms all vanish except those in u^ 9 v^ 9 iv^ 9 % 2 , v^ 9 w^ 9 u 1 u 29 v l v 29 
w x w 3 . Also from symmetry the coefficients of v^ 9 v£ 9 v 1 v 2 , are equal respectively to 
the coefficients of w^ 9 w^ 9 iv L w 2 . 

In what has gone before we have expressed the coefficient of u^ 9 %b^ 9 u x u 2 in terms 
of quantities determined by the radii and distance of the spheres, and have shown 
how the coefficients of the other terms depend on the images of the motion, whereby 
we can without much difficulty approximate to their values when the distance of 
the spheres is large compared with their radii — or the distance between their 
surfaces is large compared with the radius of one of the spheres. We pass on to 
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consider the case where one sphere is inside the other. An approximation to the 
value of the coefficients of v-f and v l v. 2 is given in § 15. It is remarkable that in the 
case of two cylinders the coefficients of the terms in v? 9 v* are equal, while those of 
u Y u % and v l v 2 are equal and opposite. But this is due to the fact that in a cylinder 
the image of a doublet (or a source) is a single doublet, whatever be the direction of 
the axis of the original doublet. 

9. If S S x be in a line through the centre the infinite trail of images of § 2 cuts 
out, and we are left with an image source and sink, and a line sink between them, 
supposing S to be outside S r Let now S and S x approach together and become a 
doublet whose strength is /x. Then we shall get a single doublet as its image whose 

Lbd S'S' ux$ 

strength =Y-L.-^ i = — — as in the former case. This we might have deduced at 



bbi 



¥ 



once from the case of the external doublet in § 4, considered as the image of its 



image. 



Fig. 5. 




If we proceed to find the kinetic energy, as in the previous case, we must clearly be 
led to the same/orm for the result, viz. : when A is moving with a velocity u from B 



2T=p£ 1 tt a « 



1 + 321 



00 / Ml 



fa 



p, 



where \L n . . . . are the strengths of the doublets inside A alone. But in this case the 
relations between the /x, p, cr are given by the equations (a, b being the radii of sphere) 



l^)i — 



v 



>n- 



n' 



a" 



cr, 



whence 



a° 



b 



3 



C+(T n = 






C+ p 



n-\ 



a" 



Pn— £ 2 



C + pn-i 



— C 



C?-W 



a 



2 



p 7 ,pn-\ + " pn + — -p«-i + CI 



2 =0 
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which differs from the equation for external spheres in having — p for p for all values 
of n. We may therefore use the same solution and writing here 



OA= v / X 8 + 



a 



2 



v / X 2 +6 3 -^H^=c 






^_ & 3 + c 3 — a? 

QB= — =7-0 



c=r — - r 



9 1 



a 



a 



r x + \ 



r z —\ b 



b r % + ^ 



q= 



Si 



p«=(^i-X)j 



1-5 



2'M 



^Y« 



which is the same form as before, only q is the inverse of its former value. 
And, as before, 



2t=|m 1 # l+sa-ft^ru- 



=iM 1 ^{l+3Qfe.ft)} 

A table for Q when 6= 2a is given at the end of the paper. 

10. It will be well here, before passing on to the consideration of the motion, to 
make a short digression on the properties of the functions Q and QV In the first place 
it is easily seen that the series for the Q and Q' functions are both convergent, even 
up to the case when the spheres touch, or q=l ; for the ratio of the n th term to the 
n — 1 th is 

r wry* J 

and this is always less than q 3 , which is less than unity, except in the case when the 
spheres touch. In this particular case the n ih term tends to the limit 



1—n 



dq 
'dq v 

3 P 2 
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do 

and the series is still convergent. The value of -r~ is the limit when \=0 of 

< h~ ( h _._ a + i 



Hence when the spheres are in contact 



q=s: 



n(a + b) -f b 



|3 — ^,3 






1 \3 



n 4- x. 



= -^X5 1 °g r ( 1 '+»). . . . (9) 



if x= 



b 



The values of this may be found from Legendre's table of the 



log r-functions. 

If the spheres be equal x=^ and 



Now 



whence 



1 



vJ — — 2Lf i ,^ -i \'i — ^ S *"""" * * * * * * * * \ / 



(2n + lf 



1 . 1 



S 3 =l+-+-+ .... =1-202056903159 . . . 

Li O 



S / 3 =|S 3 = 1*051799790264 



i . i 

When the spheres are equal g r 2 == ""* ^ in ^ s case ( 1 denote either q^ or — 

q 1 ■ q Y 



Q=(l-2 s ) 8 2r 



2 



2m \ 3 



l-04»+3 



11. We may easily express the general term in terms of r, a, b. 
in the form 



For writing it 



Q=< 



r i 

2i 

J] 



ll _ f 

2* ?i 



>=u/ suppose, 



we get at once from the relations (2), (3), (4) 



27ut u b 



a. 



IU= 



j u 



(r 2 + X) w+1 (r 2 + X) w - (r x — X) tt+1 (^3 - x ) 



which, since r%=r — 1\ and r-^— X 2 =a : 



2Xa>% u 



( r i + x ) M r i + ^) "~ ft3 } w — ( r i ™* x ) M r i — ^) 



■a a }» 
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No^ 



rr l -—af=- 



r 2 + a 2 — b 2 



— a 2 = 



r 2 — -a 2 — ft 2 a? 3 



9 =7^ suppose, 



also 



W = ^2 _ a 2^,2 _ 



(r 2 + a g -5 2 ) 2 -4ft 2 r 8 _^-4a 2 & 3 

_ 



We shall further write 4a 2 6 2 =a 4 
Then 



2\ci n b n 



u u — 



r x i 



x z 



• rffi \n 



X A 



+ \r )* — ( ^ —\r) \ +\\ (-- +\r )* + ( - — Xr 



9, 



/ i l\ 2 

2 u a n b 



x 
2 



2 \H 



2rr^t< 



2j9 + l w— 2p — 1 v 7 



+Si 



n 



2p k — 2p 



x ^^ P (2\rfp 



2 n a u b 



n 



2 n a n b 



it 



arid 



Denote 



by S„.0 
Then 

Let 



Then 



71 



2p + l n — 2p 



{ (n+l)x*+2(n—2p)ci 2 }x %l ~ i P~\x*—*y' Vn 



. . (11) 



v«=2S 



w b (— i)' 7 



bind in general 



2p + l n — 2p \q \p—q 



{n+ lx 2 +2()i-~2p)a 2 }x^-^-Za^ 



n 



P 



ij)=q 



2p + l n — 2p — Yq p — q 



V „=S-I|nS« +li? .cc 2 «-^+2S«. f/ a 2 ^-^- a }a^ 



j/ _ (l + V / ^-(l-v / ^ _ = ^ | n{n-l)(n-2) x | 



2^/^ 



+ 



72, 



2p + l w — 2p — 1 



x^+ . 



S w#0 = value of 7/ when a? is l = 2 / '~" 1 
S ;a = value of ~f~ when a? is l=-(w — 1)2 M "" S 



& 



S„«= value of 



1 <% 



when 05= 1 



Now g<n. Hence 



— l__^l_z = o when ,x= 1 



clx r i \ /, 



V 



# 
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Therefore 



S 



ii,q' 



1_ 

'2\q 



_dx2 ^/x 



^=1 



Also we see that v n is a rational integral function of (r 2 —cft—b 2 ) 
When one sphere is inside the other the series for Q is still convergent. 
12. When the spheres are concentric 



Ty=.r^==. oo X= oo 

q, a 



ar 



id 



Q= 



a 



3 



Of 



\ — 



'a\s— j3_0S 



whence 



^ # d — &° 



9 • » 6 S » 



. . . (12) 



which agrees with the result found by Stokes in his paper of 1843 before referred to. 
When the inner touches the outer \=0 and 



eP 



Q — x s % . x 3 — — ^a; 3 „ log,, r(l +as) 



(n + x) 



dx d 



= -ri5129^ 3 — log 10 r(l+a?) 



^ 3 



• 6 



where 



. • (13) 



cc= 



&— & 



If as is an integer = m say 



Q=m 8 ]S 8 -S 1 






y =m 



3 J 



2020569 



— % m ' 1 [ 



a finite expression, and in this case 



m — 1 7 
a=- b 



In the particular cases 



a=^& 



Q= -61645 



a=f6 



Q = 1'08054 
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If x is of the form 



2m +1 



Q=(2m+1) 3 



O £4 



3* 



ml 



\2n + l 



= (2m+l)M -0517998— VI 



2n + 



i) s l 



Also a finite expression, and in this case 



In the particular cases 



2m — 1, 

a= A 

2m + 1 



a=£6 Q= '39859 
a==f& Q= # 84535 



The expressions for Q directly in terms of r 9 a, b are the same functions of 
a 2 +?> 3 — r 3 as the corresponding expressions for external spheres are of r % — a 2 — b 2 . 
13, The series for Q' may, as in the case of Q, be shown to be convergent 
When the spheres are in contact 






ah \ 3 



3 



a-\-h. 



k5q • 



• * 



(14) 



Also the general term in r, a, 6 is given by 



^ 



(2a&) w 



2 n -~ 1 a n h n 



w 



n 
2£>+l n— 2p — 1 V / 



rS( — l)^^^^- 4 *" 2 



• • l JL %) J 



It is easily seen that both Q, Q / for external spheres diminish as x — i.e., r — 
increases. 

dQ dQf 
Hence for external spheres -7- , -7-7 are both negative. 

1 dr dr ° 

When one sphere is inside the other, Q decreases as x increases— i.e., as r diminishes, 

Hence in this case -7- is positive. 

ar x 

The values of the first three terms of Q, Q' are 



forQ 



and for Q' 



ah I s 



/ItA _ _ Jyi 



w 



% 3 



3 



( / y" 11 n \ — — fi™ / y > *' 



873 



a s & 



a? 6 + a 8 a5* - 2a% 2 - a%% 



'ab\3 /a}W\^ 



a?¥ 



v. 



• • 



r ) \ rx? J * \r(x 4i — a% 2 ) ^ 



3 



J 



. (16) 
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14. We may easily find jf, — at contact of the spheres. If Q* denote the w th 
term of Q x , then it may be shown that, x denoting 



a + b 



dQ n __ n(n+l)(n — l + 3x) 3 

dr a{n-\~xY ■ 



c£Q' w % 2 — 1 + 3#( 1 — #) o 

__L____ i- /'yj.v 

— g ^ 



aPdr an 

both of which are of the order -. Hence the values of -r 1 , -p, — at contact are 

n dr dr dr 

d 1 

= — oo. But though this is the case, the value of — (Q a — 3 Q') at contact is finite. 

The n th term is 

n* (n + l)(n- 4) -f 3 (1 -#) fw + a?) 4 + O 3 - 1)(6# + 4m + ar> 3 

- . ^__ — .._ . . -. - - ~- -- - -- — - j^ 

1 . _ 

which is of the order -=, and therefore the whole sum is finite. Also when n>2 the 

7i th term is positive, when n=l the sign depends on the value of x. But by consider- 
ing the values of Q, &c, in terms of r, expanding them in ascending inverse powers of 

rJ / 1 \ O' 

r, it can be shown that ~i Q—^Q I is positive always. Further, at contact Q is 

dr\ a 6 ] x co° 

a negative quantity, whilst at an infinite distance it is zero. Hence, on the whole, it 

7 / 1 \ 

must increase with r, and if this takes place continuously, y(Q — ^Q j wou ld always 

be positive. Though I have convinced myself that such is the case, I have not been 

able to prove it in general. When the spheres are at a great distance the values of 

1 
Q and Q' depend only on their first terms, and Q—^Q' only on the term of Q', which 

is of the order -. Hence here also the differential coefficient is positive. I have cal- 

culated and laid down curves representing the magnitudes of the Q and Q'-functions 
in the case of equal spheres, and when the radius of one sphere is twice that of the 

7 / 1 \ 

other, and in both cases the value for T.IQ - lQ') comes out positive for all distances. 

In what follows we shall suppose that this quantity is always positive, but it must be 
understood throughout as on]y proved, for the case of equal spheres and the case in 
which the radius of one sphere is double that of the other. 

15. Although the rapidly increasing complexity of the successive images when the 
spheres move perpendicularly to their line of centres would lead us to regard the 
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problem of finding the energy in this case as almost hopeless, yet we can carry the 
approximation to any number of images with less labour than might at first sight 
appear. For suppose we wish to take into account 2n images in A, due to A's 
motion, that is on the whole 4?i reflections. We need only first calculate the distribu- 
tion of doublets for a general position of the original one, in the ?i th image in A, and 
find the amount of the first n images. We can then treat the second portion of the 
2n images as the images resulting from the different parts of the n th image, and 
employ our first result to find the amount of the second portion by a single integra- 
tion. Suppose we proceed as if we were going on indefinitely : we suppose an original 
doublet in A at a distance p and calculate the density of the parts of the first image 
in A, say/(r) at a distance r, and thence its amount. We employ this result to find 
the density at any point of the second image, regarding it as made up of images of the 
different parts of the first, and this we do by using the expression found before, sub- 
stituting for the original doublet at p, an amount f(r)dr at a distance r, and integra- 
ting with respect to r over the first image. Thus we find the distribution for the 
second image and its amount, and therefore the amount for the first two images 
together. Starting now from this, and proceeding in the same way, we find the dis- 
tribution and amount of the first four images, then of the first eight, and so om 

Thus to find the distribution of the 2^th image we only require jp+ 1 operations, and to 
find its amount only p operations. Even with this method of proceeding the work 
would be exceedingly laborious. But for all practical purposes the first two images in 
A, i.e., the motion due to four reflections, will be sufficient— except when the spheres 
are in contact. We proceed then to find the values of the coefficient of v^ and of v y % 
to this degree of approximation. 

Suppose we have at P inside A a doublet k at a distance p Y from A, whose axis is 
perpendicular to A B. 

i. First image in B. — Then we have at Q l5 its inverse point in B, a doublet (^7) h 

k 7* 

and a line doublet thence to B, whose line density = — T'z^zr. 

J b BP X 

ii. First image in A. — The image of this in A consists of two parts, that depending 
on the single doublet in B, and that depending on the line doublet. 

(a) Image ofQ h — A doublet at P 2 ( AP a = -r^r- J whose magnitude is ( t~-^: ) h, and 

k I b \ 3 R 

a negative line doublet from P 2 to A whose line density = 'Pop") 77T 

(/3) Image of negative line doublet.— -At a distance r from B we have a negative 

k r * . ($ 
doublet = — t'^t <^ r - This has (1) a negative doublet at a distance from A= =R 

Br \ v ' ° c— r 

(ct \ k r 
-— ; ) T'^f dr. That is from P 2 to A we have a negative line doublet 

whose density at a distance R is 
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a 



k 



O //' ni r/rf 



ell 



fi~~rl bBF^lll 



and c—t 



R 



dr 



a" 



<£R"~R 3 



.'. density 



k cR—a? 



(2) a line doublet image of each portion. The doublet 



a* 



k 7 



b Bl\ 



~ dr produces a positive 



k rdr R 



line doublet from a distance R'=- — to A, whose line density ~ _ __ 

c-~~ r J ab Bl\ c—r 

Hence the density at a distance R, due to this part, from the whole line doublet 
in B 



ibq; 



k R rdr 



E 

k R 



a*ab BP X c — r 



ab BP X 



clog 



R \R 



1_ 



Hence finally the density at a distance R from A of the resultant line doublet 



h I b \s R 



A cR—a 3 , A R 



a XEVJ AQ X ab BP 2 ' aft BP 



logf- 



,21 



1 



■on- — , 



h B 



& \3 



a AQi VBPJ ' aJ BI 



+3:W"iclog 



ft 

E 



BQ, 



bind the whole amount 



6' 



* f p * 1 4, 0-Y +S- — w s ^ Rem 



ft 



aJoLAQjVBlV ' &.BP! kBl-^ 



<x& 



BQj.a 8 



E 



a 3 c 



^yAQ^BBj/ " 1 " S 5.BP 1 .AQ 1 8 " *6.BP 1 .AQ 1 9 . 



So that the whole amount of the image is 



*• \AQ r BP J k+ H.m\ AQ^ C ™ 2BQl ^ /c 



or substituting for AQ 1} &c, in terms of p 1 



*\(a-&- e pJ *" t "46(c»-6 s -cp 1 ) 8 
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For the first image in A, pi=0 
Hence 

PI — ¥] 



ab 



|3 



a?{c%-2W) 



c % —ffi J l *>7»/^2_.7»m 



2b(e*-bJ 



IV m 



• •••*** 



(17) 



The density at any point of the first image is 



Tcb 



1 , p% 



a M-S* 5 s fe E 



log 



kAkV 



8 — p J A; at a distance p % 
The amount of that part of the second image in A which depends on the latter is 



together with a doublet ( ^^ ) h at a distance p. 



i1 y _j»_ Cft J +2&(^&^^ 



&& \ 3 



c 2—1)% 



i 



and the amount of the part due to the portion of the former at a distance R is 



hh\( ab \8 , a3( c s_ 2 j2_ c R) 






whence the whole amount due to the former 



a%{c^W) 






aV 






^JoLC^-^-cB) 8 2(c 8 -& a -cK)» ' 2&(c»-&»-dft) 26 



c »_j8 &»~&R 



log *~ [ <$R 



ka?b 



V{$-~W) 



4^(cP-6*)L(<*-J»-cR)» <*-6*-cR b b 






P» 










&0- 3 



lOff v 



,ft 



lb 



46* Jo c*-8*-cIfc 



rJ"R 



/ca 3 p»., p, d j" 5 8 (e 2 — 5 2 ) 
+ 4^1 l0 ^ " "" 



5(c 2 +5 2 ) 



Hr-^dR 



E ^EL(c 2 -& 2 -cE) 2 c 2 -5 2 -cE & V-J 8 



fc 2 5 4 



| 



Z?a 2 



4 C 3( C 2_ & 3_ C/J ^2 I 4( fi 2_J3) 



2& 2 



25 2 -c 2 /;> 2 (5 2 -2c 2 ) 



. rt1 c 2 — 5 2 — cp 2 l fe 2 /* 1 log* , 

+2 log — -■■ " >■ — —I -x— v^ — cfe 



e 2 — 6 s 



463 * *=» 

% 






W 



Now ®p%—~YZ3* = <%v a Bay. 
Then the above is 



3 Q 2 
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W 






+ 



ha? 



4c*\l-«V ' 4c 3 



2& s 



1-a 3 



+ 2 (a 3 



.^+^+2(c»-6»)log(l-«») 



&a 2 r* 1 log x , 
46»\ 1 



/fa. 







2 \3 



S\ 1 _ r .3 



'^ 



iff+ ? 



/3 



2 



$ \l-~a? 






+ 2 



y8* 



1 



a 



2 



, , / a? , 6 2 \ O0 , ^ , <x 3 a 3 1 ,., ox a 3 P 1 log x 7 



45 3 



a 



2 



•a? 



6a 



Wherefore the whole amount of the second image in A x is (writing ^==-^— 



b* 



. (18) 



Ms 



Mi 



Substituting for — and ™ in (7) we get the part of T depending on v^ correct to 

the second image in A. Interchanging a and 6, the part depending on v % 2 is found. 
In the case of equal spheres 



Mi_x 

Mo 2 



#' 



fl—a? 



3 0g(c8--2& 8 ) 

+ 2(c*-a 8 ) 8 J 



M2_ i_P3_ jl ft4g2 +( c ' a / pa_Li 



a' 



Mo 
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a' 
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4' 






r4 



4\f,2 ,,2N2 



(c 3 -<)' 



+ 2P 



<r 



(c 3 -a 3 ) 3 c 3 ' 2(c 3 -a 3 ) 



^2 



log(l 




lk ^cfa 



' t-v 



16. To find the value of the coefficient of the term in v l v 2 we need to find the 
amounts of the images in B due to the motion of A, and vice versa. 

/ b \ 3 
The first image in B of k in A at a distance p l is f — - J h at Q l and a negative line 



k 



BP L 



doublet thence to B, whose line density is — "r-pp"- 
The whole amount is therefore 



6 \ 3 / b V 



fc BP x 



2 



BP, 



2 



* 



Pl 



^ /J\8 

For the first p!=0 and ~ I =i( -) A 

To find the amount of the second image in B we start from the first image in A 

(ctb \ 

doublet thence to A, whose line density is 



b.BF x \ S E e 



1 / 5 \s 



&R 



^ 
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The amount in B from the former is 



l 

2 



affi 



l(c-ft)(c»-P) 



h 



and from the second is 



So also 



b \» 



9 JoV«--B, 



1, /> 3 5 



log«- 



R c 3 -& 3 



*B 



a 



<m 



4ac 1 c 



ft 



Pft 



BS| ft j_n n ,/~ft__ "^ 



Aft 9 
4ac 

**! 

4$# 



Cft c 3 -& 3 -c/) 3 



_c 2 — 6 3 ( c — ft)' 



H l °g 



c --ft 



a 3 (e 3 -~ & 3 ) 3 — a 3 o 3 j , 



_ c 2_ & 2 ^_^2_ & 3)S 



1- 



a 



9 



c 



3 ~2> 3 






a6 a 



c(e 3 —- & 3 —- 6 3 ) 



3 J2 

+ uS • • • • ) 



— 2 



I>, 







a 3 & 



^(c 3 — $ 3 ~- & 3 ) 



3 i a ' 
-J-. — •< 



4&c 



_ c 3 - a 3 ' ( c 2-a 3 -6 9 )3" "*" ^ g ~"^ 



c*—a* 



Whence from (8) 



L'=2^(«^) + w(^ 



■#M' 



'a&\ 3 

7/ + 



a 3 & 3 



otLm 



_c(c 3 — a 3 — & 3 ) 3 
a 3 Z> 2 



3 a 8 & 3 r (c 3 -~a 3 --5 3 ) 3 

8c g (c 3 -a 3 )(c 3 -& 3 ) 



£4 



c 2 (a 3 +5 s ) ~2a^ 2 -~c 2 -^"1+ 



# 



4 



c 3 — a 3 c 3 — -& 3 , 



. . (19) 



Similarly can be found the coefficient of v 2 when one sphere moves inside another. 



Motion in the line of centres. 



1 7. When the two spheres are moving in the line of centres the kinetic energy is 
given by 

2T=A 1 ^ 1 3 +A 3 % 3 — 2B%t% 
where 



A 1 =m 1 +im / 1 jl + 3Q^.gj| 



g =r 3ivr / rVi 



M'Q'te) 



482 
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and m x , m/ 5 M' respectively denote the t mass of the sphere (A), the mass of fluid 
displaced by it, and the mass of fluid in a unit sphere. 

It is to be remarked that A l9 A 2> B are functions only of the distance between the 



d . cl 



spheres, and that therefore j~+t~ =0. Since no forces are supposed to act on the 



ttw-j tVwct 



system, both the energy and momentum are constant. Hence 



2T= constant 



« +« constant =d > 

0%-. OUct 



(20) 



-/ 



bT bT 

The last equation also follows at once from Lagrange's equation since -7+77= 0, 

Oil 1 Othiy 

and may be written 

(A } — B) tij + ( A 3 — B) u % = d 

We shall transform these equations by referring the motion to the velocity of an, 
arbitrarily chosen point P between the spheres, and the distance between them. 

Let P divide the distance (r) in the constant ratio -— =—. Then if x is the 

x 7 1— -a p 

distance of P from the origin, u its velocity 



C4J..1.UL 



whence 



(A x +A 



2 



OCl «— — • QC "~p" 0t/ j eA^O — — ■ £v **""" £* / 



a • 

Wjznw+ar, u 2 =u — fir 



2B)m 3 + (A ia »+ A a ^ s +2a/8B)r s 

+ 2{a(A 1 -B)- / 8(A 3 -B)}W-=2T ^. 
(Ai+Ag— 2B)w+{a(Aj— B)— ^(A a — B)}r=cZ 



(21) 



which we shall write 



whence 



or 



pu 2 +qr%-\- 2lur=. 2T 

* 

pu+lv=-d 

(pq — P)r^=2Tp—d % 






a! 3 
IP 



* * » * * 



« * f ZiJi) 



in which we are to take the positive or negative sign according as the spheres are 
separating or approaching one another. The spheres will move as if they repel or 
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attract one another relatively according as y « -r—f — ™ r * s P os iti ye or negative. This 

condition does not depend on their relative motion at any time, but only on their 
distance and the ratio of the constant energy to the constant momentum. The above 

condition may also be expressed, writing ~~=P, as the sign of 

P|(A 1 A 3 -B 3 )-{(A 3 -B)^+(A 1 -B)^+2(A 1 -B)(A 3 -B)f 

The last term is positive, for A 1? A 2 , B all decrease as r increases. Now h must 

always be <p since r is always real. If we put F=p=A 1 +A 3 — 2B in the above, 
the criterion reduces to the sign of 

(A 1 A 3 -B^)|{A 1 +A a -2B} 
i.e., since A x A a — B 2 is always positive to the sign of 

|(A 1 »B)+|(A 3 -B) 

d 
Now we are led to conclude from the argument in § 14 that —(A l —B) . . . are 

cor 

always positive. Hence when k has its greatest possible value the criterion is positive, 
much more then is it so for any other value of k Hence we are led to conclude that 
whatever be the relation between the momentum and energy the spheres always move 

so that r tends to decrease, whilst in the case of equal spheres, or that in which the 
radius of one is twice that of the other, we know for certain that such is the case* 
We cannot prove from this that the spheres move with reference to a fixed point as if 
they repel one another, for it might happen that both the spheres might be accelerated, 
the extra energy of the motion of the spheres themselves being taken from the fluid 
motion ; or that both are even retarded. We can easily show, however, that both 

cannot be accelerated if r is positive and both move in the same direction, for the dis- 
tance in this case increases, and therefore so do A l — B, A 3 — B, and hence because 

(Aj— B)%+(A a — B)% is constant u l9 u % cannot both increase. Also if r is negative 
and u l9 u^ of the same sign the same result holds. 

In the case where the spheres are projected. so that the momentum is zero 



• 2 2Tp 



A^-B 3 

and the relation between the velocities of projection that this may be the case is 
given by 



484 MR. W. M. HICKS ON THE MOTION OF TWO SPHERES IN A FLUID. 

U-i X3L{> ~~~ Jt> 

lOq An ""~~ JD 

When the spheres are equal %=•— % and the motion is the same as that of a single 
sphere in a fluid bounded by a plane, and moving perpendicularly to the plane. 

For this particular case 

4T 



ni«— - 



A + B 



or if to denote the velocity relative to the fixed plane r=2i^ and 

o t V^-"^"-^/0 9 

U ~A + B~ A + B U ° 

where (A+B) , u are the values of A+B, and u at the point of projection. If the 
sphere is projected from contact with the plane 

(A+BJosm+^m'+fm'diSg— l+-g-S s ) 
=m+£m'+ -3030853m' 

= m+ "803085m' 
and at an infinite distance 

A+B=m+^m' 

Hence the ratio of the limiting velocity to the initial velocity is 



v 



1 + -6061707 l 



2p + l 



where p is the density of the sphere. 

For densities 0, 1, 10, the values of this ratio are respectively 1*2661, 1*0963, 
1*0143. The greatest value is when the density of the sphere is zero, and the least 
is when m'=0 (no fluid) or m= oo , the ratio then being, as it ought to be, unity. 

In the case where the spheres are unequal and projected with no momentum from 
contact their initial velocities must be opposite and in the ratio of the quantities 

m s +|m^-|m^ 3 {lD 3 log, r(l +*/) + S 3 } 
and 

m i + i m 'i — fw'ifc 3 { |D 3 log, r( 1 + ®) + S 3 } 

h & 
x and y denoting the quantities — — , — r 

CO ~t~ CO "j" 
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If a=2h, #=■§-, y=f, and we find from Legendre's tables "of the Eulerian integrals 

D 3 log 10 r(l+^) = -'485 

D 3 log 10 r(l+2/)=--275 

and the ratio is 

1 P + -1174 

8 p + '4642 

which when the densities of the spheres and fluid are equal becomes 

* 763 Ann, 

-— = •0954 

We find the velocities of the spheres relatively to the fluid by eliminating u 

between 

» 

and 

pu-{-lr=d 
whence 

__rf A 3 — B ^ 

p p 

and 

<$ A,— B * 
w 9 = — r 

Suppose now the same spheres projected with the same initial circumstances except 
that now the spheres have changed places, and let u 2 , u\ be the corresponding velo- 
cities at the same distances. Then 

, d , A n — B • 

a 9 = — |— ^ r 

p p 

since d and r do not depend on the question which of the two is foremost. 

Now if a > b we see at once from the expressions given for A l5 A 2 in terms of the 
distances that A x > A 3 , and hence that the foremost will be most accelerated when it 
is the smallest. 

If now u l9 u% denote the velocities at any moment which we may regard as the 
velocity of projection 

2T A l tt l 2 + A 2 itc^ — 2Bu l u 9i 

Writing £ for the ratio - 1 the equation to find £ in order that k may have a given 

% 
value, is 

MDCCCLXXX. 3 R 
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This enables us to find within what limits k must lie, for f must have real roots, 
and therefore 

or 

Hence W may be any positive quantity less than p. The greatest possible value of 
this is when the spheres are infinitely distant, and then 

^=m 1 +m a +-^(m' 1 + m 'a) 

To each value of <f will correspond two states of motion, the initial velocities in each 
case being opposite. For example, if f is positive, ie,, both velocities in the same 
direction, the two states will be when (a) is the foremost, and when (6) is the foremost ; 
if ^ be negative, the two states will be, one in which the balls begin to move towards 
each other, the other in which they begin to move from each other. Thus for every 
given value of k there are four possible states of motion. 

If ever ^=0 then ^=0, and the spheres must be at such a distance that 

Now, supposing k given, this can only happen if k % lies between the greatest and 

(A. --B) 3 
least values of ■ 2 . ■-. The least value is when the spheres are in contact, the 

greatest when they are at an infinite distance, the value then being ,w 3 +2 m V 

. (A, — B) 3 
If ^ 3 =0, then F must lie between the greatest and least values of ~™~~— . Now 

Ag N A x 

as 

(A^— B 2 )(A 2 — A^O 

as 

If we suppose a > b then A 1 > A 3 , and calling kf 9 k^ the least values of the above 
limits &i<&. 
Hence it 

Kf<h or >m 1 +^m\ 



v 
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tlie spheres can neither ever come to rest ; if 



Vd 2 



the small sphere can never come to rest. 



The effect of the fluid on vibratory motions. 

18. Suppose each of the two spheres attracted to a fixed centre of force where the 
force varies as the distance. Let x v x % be the distances of the spheres at any time 
from their respective centres of force measured in the same direction. Then 

2T = A-^-j 2 + A 3 w/ — 2 ~Bu x u 2 = C — ni^x^ — m^x^ 

Also since we neglect squares of small quantities in finding the small vibrations, the 
equations of motion become 

» a » • 

Jl\.tX] ~~~ XjX& -— — ■" iii^LLjOu-i 
• • • • 

mmmm XjX-\ ~~\ **■■*- 2 2 """"" """""" ^^qI^QiQ 

and we suppose the spheres so distant, and their motions so small, that we may neglect 

the small changes in A, B during the motion. The spheres must not be too close, for 

dA 

at contact -~, &c, are infinite, as was shown in 5 14. 

dr ° 

Solving the above equations in the usual manner we find 

£^=1^ sin (K^+^+Nj sin (K 3 £+/3) 
x 9 =e'L l sin (K^+^+e'Nj sin (K 2 <£+/3) 



where 



Kx 3 



2(A 1 A 3 ~B 2 ) 



_ A 1 K 1 »-m 1 j* 1 _ BK^ / _ A 1 K 2 2 -m lA 6 1 _ BK 9 8 

— BKi 3 — AjjK^-m^ — BK 3 3 ~A 3 K 3 3 --m 3 /x 3 

From this we see that, to the first order of small quantities, the mean position of 
the spheres is not altered, or to that degree of approximation there is no mean 
attraction or repulsion. 

If we regard the spheres as two pendulums swinging in the fluid, in the same 
horizontal line, of lengths l l9 l 2 , then the motion is given by the above equations if 
we write 
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„ — Pi- 1 1. „ ~Pizl.lL 

where p l9 p 2 are the densities of the spheres compared to the fluid. 

If in the above we make m 3 = co we get the case of a forced vibration of period —7 
In this case 

1= K*=k K 3 *=-p 

o? 1 =L sin (\/fat-\-a) 

x =~ — — — - L sin (\/A^+ a )+N sin f a/^~^^+/3 



/*i 



^ 






If the sphere (b) is set free when (a) is for the moment at rest, and the time be 
reckoned from this moment 






eLl COS 4/ fat— cos A/ —r^t 



and the motion of (6) consists of two periodic terms whose amplitude is e times that 
of (a). 

Let now the strength of the centre of force on (6) diminish indefinitely. Then 

^=-r-L(cos \/fat~—l) 

and (6) would oscillate in the same period as (a), without being attracted or repelled 
towards it except by forces depending on the square of the amplitude of (a). To find, 
then, whether the action of (a) on (b) is attractive or repulsive we must take account 
of quantities of the second order of small quantities. 
The full equation of motion of (b) is 



For a first approximation we have 



B 



iX/i 



Write 



2 =Y-L{cos^fat-l) 

A 2 



B _ * 

f T 2 =— L(cos y/^— l) + « 
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where % is of the order L 3 at least. Substituting for % and neglecting cubes and 
higher powers of L, 

A(tJ\.n i JDJuyU'i i . CvXj j -j. , — 

3 • 



B 3 B \ dA __ d 

-i L Vi(x72~ 2 A )^ sin2 v / ^- + |LVi^;(A 1 -2B) sin 3 ^£=0 



and 



B\ ,_ . BL 






3/ ^ 



Whence the equation takes the form 



where 



» « 
•A a 2=/+<j( cos x/fat+h cos 2^/J^t 

2f / B\/B ^A 2 dB\ jB/B VAj i^/a o"R\ 

^Vi \ Vte^ 7 ""^/ + %\A 2 ^ 2 /^r~V Al "" 2 ^ 

i ^/ AxAg—B 

~~ 2 ^r\ A 3 

in which last form we may neglect in A l the m 1 +im\ as it disappears in the differen- 
tiation. Hence the mean action of (a) on (h) is an acceleration towards (a) 

— __ k 2 A*i dVA^— B s 

4A 3 dr\ A 3 



-£-»-q Q// \ "^-Q / 



if v is the " velocity of mean square " of . (a). 

If the distance of the spheres is so large that we may neglect twelfth and higher 
inverse powers of the distance, we need only consider the first images or the first 
terms in A and B. In this case it will be found that the acceleration to (a) is 



18u 3 /ay 



2p + l\r 



(r»-& 2 )* 2p + lr m 



(23) 



To find when there is repulsion 



i s 3 



< 



(f-b*)± 2p + i 
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or 

h 



r 



vi- </m 



which can clearly only happen if 2p+l<S or the density of the sphere less than the 

In general, then, when the body is denser than the fluid it is attracted. If its 
density is less than the fluid there will be a critical point (as mentioned by Sir W. 
Thomson), beyond which there will be repulsion, and within which it is attractive. 
This critical distance is given by 






VI'- v^ 



(24) 



in using which it must be remembered that if r comes out nearly equal to b } the 
formula fails to give a correct value, as it was obtained on the supposition that the 
distances were large. It is, however, extremely accurate if we remember that it is 
true up to inverse powers of the twelfth at least. If the density of the sphere be *9 
the critical distance would be 7 '648 times its radius. It may be noticed that while 
the principal term, in the acceleration depends on r^ 7 , if the density be the same as 
the fluid it depends on r"™ 9 . 

In the case of a sphere vibrating within another sphere, along the line of centres, 
the effect of the fluid will be represented by supposing the inertia of the sphere 
increased by a mass 

= i{l + 3Q(g'.g r ] ) X mass of fluid displaced by it 

where Q has the value given in § 9 : provided it is not close to the boundary of the 

clQ 
containing sphere, as in that case — becomes infinite, and the small motions of the 

sphere will produce great changes in the value of Q. When its mean position is the 

dQ 
centre, -^ =0 and Q may be considered constant when we neglect in our equations of 

motion cubes of small quantities. The value of Q in this case is, as has been already 

menxioneci, 

b s 4- 2a? 



JL»- 

2 ty—a 



~X mass of fluid displaced 



The foregoing serves to solve the problem of a ball pendulum within a spherical 
envelope when it is so suspended that its centre lies in the horizontal line through the 
centre of the envelope. When it oscillates in any other position the value of the 
coefficient of inertia may be approximated to as in §§ 15, 16 t 
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19. If instead of supposing the sphere (h) free to move we suppose it held fast, and 
require to find the force necessary to do so, we get a different result from the fore- 
going. Suppose the sphere (a) moving in any manner, the sphere (6) being for the 
moment at rest, and suppose a constant force F acting on 6. 

The equation of motion for B is 

* * * * d .A. d 

A^-B^+W^i-%)^T+W^(A 1 -2B)=F 

Suppose now that F is of such a magnitude that ^6 3 being zero it makes x 2 also zero. 
Then F is the force required to keep (h) at rest at the moment when the motion of (a) 

is given by u v x y Hence 

P=-Bx 1 +|V|(A 1 -2B) 

Let x l ='L sin K.t, L being small. Then neglecting cubes of small quantities 

F=(b+^ clAhK* sin Kt+± j T (A 1 -2B)L 3 K 2 cos' 2 Kt 

and 

dr=x l ='L sin K.t 

.\ F=BLK* sin K*+iL 8 K 8 j^+i|;(A 1 -2B)| +^cos 2Kt 

This is the force at the time t necessary to keep (b) at rest. Hence the mean force 

dA . 

is a force =^L 3 K 2 ~-~ towards (<x), which is equal and opposite to the force of (a) on 

dA . • . 

(6). Since ~~ is negative, the action is an attractive one 

Oji 

~ 4 dr 

dr 

Taking for A l only the first term of Q, which is equivalent to neglecting twelfth and 
higher inverse powers of r 

ab \ 3 



A l =m l +^m\< 
and the force 

= 9mV. cm ' 



1 + 31 



/ytZ _ _ J\Z 



i)* 1 afar 
^^-•t^t^X weight of fluid displaced by a . . . * * (25) 

is \ / 

For example, for equal spheres at a distance Aa (distance between their surfaces 
2a), the mean square of velocity of (a) being the same as for oxygen at a tempe- 
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rature of 0° G, v= 1524 feet per 1", and the force =— : - X weight of fluid displaced, a 

being measured in feet. It is clear that while the force decreases indefinitely with 
the size, the "effective" force increases indefinitely. 

If (a) vibrate through a distance yg-inch, 256 times a second, and a=^-inch, the 
force is '01197 weight of water displaced = weight of 12*8 milligrammes. 

In the same manner can be found the action of (a) on (h) when a describes any 
small curve whose plane contains (5). 

Values of Q, Q' for equal spheres. 



r 

2a 


&• 


Q. . 


a s Q'- 


1 


1 


•051800 


•150257 


1-05 


-72985 


•028307 


•116749 


1-1 


•6418 


•018768 


-098312 


1-2 


•5367 


•009531 


•073754 


1-35 


•4431 


•004049 


•0511 


1-5 


-3819 


•001959 


•037142 


1-75 


•3138 


•0007023 


•023335 


2 


•2679 


•0002962 


•015631 


2-5 


•2087 


-0000723 


•008001 


3-5 


•1459 


•0000090 


•002895 


4-5 


•1125 


•0000019 


•001373 



Values of Q 1? Q 2 , Q' when a=2b for external spheres 



T 






1 


1 


a-f-b 


Qr 


Ms* 


^ Q '- 


p Q • 


1 


•0206 


•0945 


•04452 


•35616 


1-05 


•01228 


•04298 


•03393 


•27144 


1-1 


•00862 


•02572 


•02884 


•23072 


1-25 


•003653 


•00886 


•01918 


•15344 


1-5 


•000719 


•00119 


•01090 


•08720 


2 


•000186 


•00024 


•00455 


•03643 


3 


• * 


t a 


•000016 


•00013 



Values of Q when 6= 2a for an internal sphere. 



r 
a 




Q. 

142870 


•25 


•15106 


•5 


•18046 


•75 


•25676 


1 


•61645 



